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Abstract 


This paper describes an application of the method of moments 
to determine resonant frequencies of irregularly shaped microstrip 
patches embedded in a grounded dielectric slab. For analysis , the 
microstrip patch is assumed to be excited by a linearly polarized plane 
wave that is normal to the patch. The surface- current density that 
is induced on the patch because of the incident field is expressed in 
terms of subdomain functions by dividing the patch into identical 
rectangular subdomains. The amplitudes of the subdomain functions, 
as a function of frequency , are determined using the electric-field 
integral equation (EFIE) approach in conjunction with the method of 
moments. The resonant frequencies of the patch are then obtained by 
selecting the frequency at which the amplitude of the surface- current 
density is real The resonant frequencies of the equilateral triangular 
and other nonrectangular patches are computed using the present 
technique, and these frequencies are compared with measurements 
and other independent calculations. 

Introduction 

Microstrip patch antennas have been studied extensively in recent years because they offer many 
practical advantages, such as being lightweight, able to conform to the body of a host object, and 
relatively easy and inexpensive to fabricate compared with other types of antennas. As a result, 
many analytical approaches have been proposed to analyze microstrip patch antennas (ref. 1). The 
electric- field integral equation (EFIE) approach in conjunction with the method of moments is the 
most widely used technique to study microstrip patch antennas (refs. 2 to 6). Early work (refs. 2 
to 6) on the application of the EFIE approach to patches is limited to microstrip patches of regular 
shapes, such as those that are rectangular, circular, or elliptical. Furthermore, these earlier works 
(refs. 2 to 6) use an approach in which the surface-current density on the patch is expressed in terms 
of entire domain functions. This paper describes the EFIE approach in conjunction with the method 
of moments for solving the problem of irregularly shaped microstrip patch antennas by expressing 
the surface currents in terms of subdomain functions. 

Several authors (refs. 7 and 8) have used the EFIE approach to analyze irregularly shaped 
antennas. The work in references 7 and 8, however, is applicable to irregular plates in free 
space. In these papers, a nonrectangular plate is viewed as an interconnection of quadrilateral 
plates. The currents on the quadrilateral plates arc expressed in terms of nonrectangular surface 
modes. Although the segmentation techniques used in reference 7 completely fill the area of the 
nonrectangular plate, extra current modes are required to ensure continuity of surface-current density 
at the joining plates of the junctions. Furthermore, when nonrectangular surface modes are used 
to express the surface currents, the resulting matrix in the method of moments solution is of a 
symmetrical but not of a Tocplitz nature; therefore, the computational time increases. 

Mosig (ref. 9) and Michalski and Zheng (ref. 10) use a mixed potential integral equation approach 
to solve the problem of irregularly shaped microstrip patches. In these studies, the EFIE is solved 
by using numerical techniques in the spatial domain which may require special care to handle the 
singularity in the Green’s function. Martinson and Kuester (ref. 11) use generalized edge boundary 
conditions to accurately analyze irregularly shaped microstrip patches. This approach, however, is 
valid only for thin substrates (e.g., dj Aq < 0.01). 

This work describes a segmentation technique to analyze irregularly shaped microstrip patch 
antennas. We assume that an irregularly shaped patch is enclosed by a rectangle with sides equal 


to W x and W y . After dividing W x into (M + 1) sections and W y into ( N + 1) sections, the surface- 
current density over the rectangle is expressed in terms of overlapping triangular functions in the 
current flow direction and pulse functions in the orthogonal direction. A shape function in the current 
expansion function is introduced to ensure zero current outside the patch. The shape function is 
equal to 1 if the subdomain lies inside the irregularly shaped patch antenna, and it is equal to 0 if the 
subdomain lies outside the irregularly shaped patch antenna. However, for the truly irregular shaped 
patch, the subdomain that is close to the boundaries of the patch may be partly occupied by the 
patch. In such cases, the subdomain is considered to be inside the patch if the area occupied by the 
patch in that subdomain is more than 50 percent of the area of the subdomain. This process of finding 
the shape function, however, becomes tedious and time consuming for complicated geometries. 

The EFIE equation is reduced to a matrix equation that is solved using standard matrix equation 
solver subroutines when the testing functions are selected to be the same as the expansion functions 
(he., a Galerkin solution). The surface-current density is then used to determine the resonant 
frequency of the patch. One of the disadvantages of the present method is that a large number of 
subdomains are required to achieve convergence because the edge conditions for the surface-current 
distributions are not explicity expressed. However, the Toeplitz nature of the impedance matrix is 
still maintained, thus considerably reducing the matrix filling time (ref. 12). 


Symbols 

a x * a 3/ > a z 

d 

E* 

£*0 

E s{Jx) 

E s(Jy) 

E ti 

^Xi » E 
F X S {Jx) 
F X s{Jy) 
Eys{Jx) 

Ey S (Jy) 

Eio 

E20 

Fxmn{%) y) 
F X p(x, y) 
Ey mn (^ y) 
Fyp(x, y) 

fxp* 
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unit vectors along x, y, and z axes 
dielectric substrate thickness, cm 
incident electric-field vector 
intensity of incident electric field 
scattered electric-field vector caused by J x 
scattered electric-field vector caused by J y 
tangential electric-field vector 
x and y components of incident electric field 
x component of scattered electric- field caused by J x 
x component of scattered electric field caused by J y 
y component of scattered electric field caused by J x 
y component of scattered electric field caused by J y 
electric-field amplitude for perpendicular polarization, V/m 
electric- field amplitude for parallel polarization, V/m 
expansion function for x-directed current on (m, n)th subdomain 
expansion function for x-directed current on pt h subdomain 
expansion function for y-directed current on (m, n)th subdomain 
expansion function for y-directed current on pth subdomain 
resonant frequencies of various modes, GHz 
Fourier transform of E xp /(x,y) 

Fourier transform of F yp f(x,y) 
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incident magnetic field vector 

complex amplitude of x-directed current density on (m,n)th subdomain, A/unit 
area 

complex amplitude of ^-directed current density on (m,n)th subdomain, A/unit 
area 

complex amplitude of x-directed current density on pth subdomain, A/unit area 
complex amplitude of y-directed current density on pth subdomain. A/unit area 
x-directed induced surface-current density on plate, A/unit area 
y-directed induced surface-current density on plate, A/unit area 

= V-i 

propagation constant/wave number in free space 
propagation vector of plane wave 
number of subdivisions in x-direction 
(m,n)th subdomain of induced current 
number of subdivisions in y direction 
total number of x-directed subdomains on plate 
piecewise linear distribution in x-direction 
piecewise linear distribution in y-direction 
equivalent to (m,n)th x-directed subdomain 
equivalent to (m / , n ; )th x-dirccted subdomain 
total number of y-directed subdomains on plate 
pulse distribution in y-direction 
pulse distribution in x-direction 
position vector in direction of plane wave 

reaction of p r th x-directed subdomain testing function with E XI 

reaction of p 1 th y-directed subdomain testing function with E y j 

maximum dimension of plate in x-direction, cm 

maximum dimension of plate in y-direction, cm 

Cartesian coordinates of field point 

Cartesian coordinates of source point 

= m Ax 

location of patch, cm 

mutual impedance between pth and p f th x-directed subdomain currents, ohm 

mutual impedance between pth y-directed and p r th x-directed subdomain currents, 
ohm 
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mutual impedance between pth ^-directed and p r th y-directed subdomain currents, 
ohm 


Zyy mutual impedance between pth and p'th ^-directed subdomain currents, ohm 

a wedge angle, deg 

Ax = JT+T 

A y ~ 77+T 

e r relative dielectric constant of slab material 

7 ^ free-space wave impedance, ohm 

9 t , 4>i direction of angle of incident wave, deg 

Aq wavelength in free space, cm 

Abbreviation: 


TM transverse magnetic 


Theory 

General Theory 

Consider an irregularly shaped thin patch embedded in a dielectric slab and illuminated by a 
plane wave as shown in figure 1 . The incident field can be expressed as 


E i ( x , y, z) = E,q [— a x sin (fa) + a y cos (&)] exp [(jfc 0 k, ) • r] 
for perpendicular polarization, and 

p, 

H ( (x, y, z) = — [-a* sin (fa) + a y cos (fa)] exp [(jfc 0 k*) ■ r] 
' * 


( 1 ) 

( 2 ) 


for parallel polarization, where and ko are the free-space wave impedance and the wave number, 
respectively, and (#;,<&) is the direction of the angle of incident wave. In equations (1) and (2), 
Em is the incident electric-field intensity; the quantities a Xj a yy and a z are the unit vectors along 
the x, y, and z sixes, respectively; and 


k- = a x sin (9 l ) cos (<p{) + a y sin (9 t ) sin (<fo) + a z cos (#*) 
r = a x x + SyV T a Z z 

The tangential electric field in the plane of the patch when the patch is absent (ref. 12) is then 
obtained as 


E u = Eio [— a x sin (fa) + cos (fa)] exp { j [(fc 0 kj) ■ a*] x} + { j [(fcok*) • %] y} (3) 

for perpendicular polarization, and 

E ti = E 2 o [-a* cos (fa) + a y sin (fa)] exp {j [(fc 0 kj) • a x ] x} + {j [(fc 0 k z ) ■ a„] y} (4) 


for parallel polarization, where 


Eio = 


j 2 E m cos (Qj) sin (k 0 z'k z ) 
k z cos (k()dk z ) + j cos (0 t ) sin (k$dk z ) 


exp [jk(\d cos (0,)] 
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I 


r. 

f 


#20 = 


-j ZEj 0 k z cos (0j) sin ( k 0 z'k z ) 

. cos (9i) cos ( kpdk z ) + jk z sin ( kpdk z ) 


exp [j d cos (0*)] 


and k z = e r — sin 2 (0 z ). 

Let J x (x\y / ) and J y (x\y f ) be the x- and y-directed induced electric surface- current densities on 
the patch and E S (J X ) and E s (J y ) be the scattered electric fields caused by the x- and y-directed 
currents, respectively. Setting the total tangential electric field over the patch to zero yields 


a z x [E s (J x ) + E s ( J y ) + Ej] = 0 


(5a) 


Equation (5a) can be written in component form as 

(Jar) T F X s {Jy) T F x i = 0 ^ 
Fys («/r) T Fys ( Jy ) + Eyi — 0 J 


(5b) 


where the suffixes x and y are used to indicate the x and y components of the scattered and incident 
fields. To solve equations (5) for J x (x\ 2/0 and J y (x r ,y r ), the x- and y-directed currents on the patch 
arc expressed as 


M N + 1 


J * {x',y) = E E ix (?Tl, Tl) Fxrnn ? U ) 


(6a) 


m= 1 n=l 


Af+ 1 TV 

Jy (x\ y') = E *v n ) ( x/ > ( 6b ) 

m= 1 n = 1 

In deriving equations (6), the irregularly shaped patch is first enclosed by a rectangle with 
sides W x and W y , and this rectangle is then divided into M- hi and AT + 1 sections along the x 
and y directions, respectively, as shown in figure 1(b). The quantities z x (m, n) and i y (m,n) in 
equations (6) are the amplitudes of surface-current densities at the (m,n)th subdomain. The 
expansion functions F xmn (x\y f ) and Fy mn (x f ,y f ) in equations (6) are given by 


Fxmn V ') - Pm (x ) Qn (V) 


Pymn {x , y ) — Qm {x ) Pn (?/ ) 


where 


Pm {x') = 

Qn {y) - 




Xjn-X 
X 


Ax : 


1 X ~ X TI 

< Ax 

r 1 
0 

w x 


M + 1 


((% - Ax) <x'< x m ) 

( \x m < x < (x m Ax)) 

((n- 1) A y < y' < nAy) 
(Otherwise) 


Ay 


Wy 

N + 1 


In the above expressions, P n (y') is obtained by replacing m and x' in the expression for P m (x') 
by n and y' , respectively. Similarly, Q m {x') is obtained by replacing n and y' in the expression 
for Q n {y') by m and x', respectively. For simplicity, the double summation with respect to m 
and n in equations (6) can be represented by a single summation with respect to p. If P and Q are 
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the maximum numbers of the x- and ^/-directed subdomain cells, respectively, on the patch, then 
equations (6) may be rewritten as 


Jx y f ) — YJ (p) F X p (x f , y r ) 
p= 1 

P+Q 

Jy (*', y') = E iy (p) Fyp (*', 2/') 
p=p+i 


(7) 


Using the method of moments and test modes that arc identical to expansion modes, equations (5) 
yield 

JJ E xs (J x ) F xpl (x, y) dx dy + II Exs {Jy) E xp r (x, y) dx dy + II E xl F xp , (x, y) dx dy = 0 (8a) 

where 7/ = 1,2,3..... /’, and 

JJ E ys ( J x ) F yp , (x, y) dx dy + JJ E ys (J y ) F yp , (x, y) dxdy + JJ E xl F ypf (x, y ) dx dy = 0 (8b) 


where p r — P -f 1, P + 2, P + 3, . . . , P 4- Q, the surface integrals in equations (8a) and (8b) are 
carried out over the p f subdomain. Equations (8a) and (8b) can be written in the following convenient 
matrix form: 


7P1P yP,P 
" xx & xy 

yPV 1 yP#' 
** yx ^yy 


ix{p) 

L h (P) 


v x ( p'y 

- V y{p')_ 


(9) 


where p' = 1,2,3 ,(P + Q) and Z x f and Zxy are the self and mutual impedances between pth 
and p'th current basis functions. Detailed expressions for these impedances are given in reference 5. 
The elements of excitation vectors in equation (9) are given by 


Vx ( P ) = -Eio sin (<f>i) f xp/ [(fcokj • a x ) , (fcok; ■ a,,)] (10) 

v y (pO = E l0 cos (fa) f yp , [(fc 0 kj • a x ) , (fc 0 ki • aj,)] (11) 

for perpendicular polarization, and 

Vx (p) = E 2 o cos (4>i) f xpf [(fc 0 k» - Si) , (Ar 0 k.j • a,,)] (12) 

v y {p') = E 2 q sin (fa) f ypl [(fc 0 kj • a x ) , (A- 0 kj • a,,)] (13) 


for parallel polarization, where f xp > and f yp / are the Fourier transforms of F xp r and F yp ?, respectively. 

Resonant Frequency 

The current density that is excited by an incident plane wave is obtained by solving the matrix 
equation (9). The current amplitudes i x {p) and i y {p) are in general complex quantities; however, 
at resonance, i x {p) and iy{p) are real numbers. The resonant frequency of the patch may therefore 
be defined as a frequency at which the real part of i x {p) and i y {p) is maximum and the imaginary 
part is zero. The dominant and higher order resonances of an irregularly shaped patch can therefore 
be determined by finding the frequencies at which the real part of the surface-current density is 
maximum and the imaginary part of the surface-current density is zero. 

Numerical Results 

In this section, resonant frequencies of irregularly shaped patch antennas are obtained and 
compared with the measured data and the results obtained using the cavity model (refs. 13 and 14). 
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The induced currents i x (p) and i y (p) are determined using equation (9) as a function of frequency. 
The dominant and higher order resonances of the patch are determined from the knowledge of the 
variation of i x (p) and iy(p) as a function of frequency. 

Resonant Frequencies of Equilateral Triangular Patch 

As a first example, we select an equilateral triangular patch with dimensions a = 10 cm, e T = 2.32, 
and d = z f = 0.16 cm, as shown in figure 2. For resonant frequency calculations, the patch is assumed 
to be excited by a parallel polarized plane wave at normal incidence with fa — 0° or fa — 90°. 

The triangular patch is assumed to be enclosed by a rectangle with W x = 10 cm and 
Wy — 8.66 cm, as shown in figure 2. Selecting values of M and N (e.g., M = TV = 11) and defin- 
ing the (m, n)th subdomain for the x- and y-directed currents to be inside the triangular patch if 
more than 50 percent of the subdomain area lies inside the patch, the matrix equation (9) is solved 
for i x (p) and i y {p)- The current densities i x (p) and i y (p) that are normalized to the incident field at 
the center of the patch are then plotted in figures 3(a) and 3(b) for the two incidence angles (0 t = 0°, 
fa = 90°) and (0^ = 0°, fa = 0°), respectively. From figures 3(a) and 3(b), we see that the resonant 
frequencies at which the real part of the current density is maximum and the imaginary part of the 
current density crosses zero are /o = 1.210 GHz and f\ = 1.256 GHz. 

To test the dependence of these resonant frequencies on M and TV", the resonant frequencies are 
calculated as a function of M and N and are plotted in figure 4. This figure clearly shows that 
M = TV" > 19 gives stable numerical results and that the first resonance occurs at /q = 1.249 GHz 
and the second resonance occurs at f\ = 1.276 GHz. 

Microstrip patch resonances are usually associated with the cavity modes described in refer- 
ence 14. To identify the above two resonances with cavity modes, we must plot a vectorial repre- 
sentation of the surface-current density over the patch that is excited by an incident plane wave at 
resonant frequencies /q — 1.249 GHz and f\ = 1.276 GHz, as is done in figures 5(a) and 5(b). These 
figures show that the magnitude of the surface-current density is proportional to the length of the 
vector, while the direction of the current flow is indicated by the arrow direction. Upon careful 
examination of the resonances of the TMqi (transverse magnetic) and TMp) cavity modes given in 
reference 14, it is clear that the resonance at /o = 1.249 GHz corresponds to the TMqi mode, and 
that the resonance at f\ = 1.276 GHz corresponds to the TMjq cavity mode. Note that the cavity 
model discussed in references 13 and 14 predicts that the resonant frequencies of the TMio and TMqi 
modes are identical. This technique, however, gives slightly different resonant frequencies of these 
modes; this difference may be attributed to discretization of the patch. 

To study higher order resonances of the triangular patch, surface-current densities that are excited 
by plane waves for the incident angles of ( 9 $ = 0°, fa = 0°) and (0* = 0°, fa — 90°) are plotted in 
figures 6 and 7, respectively, over a wider frequency band. The resonant frequencies of higher order 
modes of the patch are determined from these figures, and they are presented in table 1 with the 
corresponding cavity modes calculated using the cavity model (ref. 14) and measured results given 
in reference 13. 

The resonant frequencies of the TM mn and TM nm modes are identical, as seen from the cavity 
model formulations. As noted earlier, the present method that is based on discretization predicts the 
resonant frequencies of these modes to be slightly different from each other. Corresponding vectorial 
representations of the surface-current densities at incident angles of (fy = 0°, fa = 0°) and (0 l = 0°, 
fa = 90°) are depicted in figures 8 and 9, respectively, for higher order resonant frequencies. 

Resonant Frequencies of Circular Patch 

As a second example, we consider a circular patch with a = 1.88 cm, as shown in figure 10. This 
patch is assumed to be excited by a parallel polarized plane wave at normal incidence with fa = 0° 
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or (j>i — 90°. The induced current density i x (p) and i y (p) as a function of frequency at the center 
of the patch is obtained after solving equation (9), and it is plotted in figures 11(a) and 11(b). 
Values of M and TV were arbitrarily selected to be M = TV = 11. These figures show that the 
frequency at which the real parts of the current densities are maximum and the imaginary parts 
of the current densities are zero is /q = 2.760 GHz. The first resonance of the patch is therefore 
at 2.760 GHz. To test the dependence of the first resonance on M and TV, /o is determined as a 
function of M and TV, and it is plotted in figure 12. Figure 12 clearly shows that M = TV > 19 
gives stable numerical results. This figure also gives the first resonant frequency of the circular 
patch obtained using the method described in reference 15. Good agreement exists between the 
two results for M = TV > 19. The vectorial representations of the surface-current density over the 
circular patch, excited by the x - and y-directed linearly polarized plane waves at the first resonant 
frequency, are shown in figures 13(a) and 13(b), respectively. Comparison of this representation 
with the cavity model representation (ref. 14) indicates that the resonant mode at a frequency of 
2.760 GHz corresponds to the TMp cavity mode. 

Higher order resonances of the circular patch are determined, and they are given in table 2. This 
table also gives higher order resonances that are calculated using the cavity model (ref. 14). Good 
agreement exists between the tw r o methods. A vectorial representation of the surface- current densities 
at higher order resonant frequencies is given in figure 14, A comparison of these representations 
with the representations obtained by the cavity model confirms that resonances at frequencies 
of 4.685 GHz, 5.855 GHz, and 6.360 GHz correspond to the TM2i, TMq 2, and TM31 cavity modes 
(ref. 14). 

Resonant Frequency of Trapezoidal Patch 

As a third example, we consider a trapezoidal patch with dimensions as shown in figure 15. 
Figures 16(a) and 16(b) present the variation of current density as a function of frequency at the 
center of the trapezoidal patch for two angles of incidence. These plots show T that the frequency at 
which the real part of the current density is maximum while the imaginary part is zero is 1.342 GHz. 
The first resonance of the patch therefore occurs at a frequency of 1.342 GHz. To test the dependence 
of the first resonance on M and TV, the first resonance frequency of the patch is calculated as a 
function of M and TV, and it is given in table 3. Table 3 clearly shows that M > 12 and TV > 6 give 
stable numerical results. Figure 17 gives the vectorial representation of the surface-current density 
on the patch at the first resonance for incident angle of (0j = 0°, — 0°). 

Conclusions 

An electric-field integral equation approach in conjunction with the method of moments has been 
used to determine the resonant frequencies of irregularly shaped microstrip patches. Numerical 
results obtained using this approach compare well with experimental results and other independent 
calculations. Discretization of an irregular patch into symmetrical rectangular subdomains, in the 
present technique, results in a symmetrical and block Toeplitz impedance matrix. The discretization 
scheme used, however, does not explicitly enforce the proper edge conditions on the surface- 
current distribution. As a result, a large number of subdomains are required to achieve numerical 
convergence. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
August 16, 1993 
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Table 1. Calculated and Measured Resonant Frequencies of 
Equilateral Triangular Patch for M = N = 19 

[Patch shown in fig. 2] 


Mode 

Resonant frequency, GHz 

Present method 

Measured results 
(ref. 12) 

Cavity model 
(ref. 13) 

TMoi 

1.249 

1.280 

1.299 

TMio 

1.276 

TMn 

2.172 

2.242 

2.252 

TM02 

2.525 

2.550 

2.599 

TM20 

2.510 

tm 12 

3.265 

3.400 

3.439 

TM21 

3.356 


Table 2. Calculated Resonant Frequencies of Higher 
Order Modes of Circular Patch 

[Patch shown in fig. 10] 


Mode 

Resonant frequency, GHz 

.... 

Present method 

Cavity model 
(ref. 13) 

TMn 

2.816 

2.818 

TM 2 i 

4.685 

4.674 

TM02 

5.855 

5.864 


Tabic 3. Resonant Frequencies of Trapezoidal Patch 
for Various Values of M and N 


[Patch shown in fig. 15] 



N 

Resonant frequency, GHz 


5 

1.370 


6 

1.345 


7 

1.342 
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1.343 

18 
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1.343 

20 

10 

1.342 
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(b) Top view of microstrip antenna. 
Figure 1. Geometry of microstrip antenna. 
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Resonant frequency, GHz 
(b) (0i = 0°, 4h = 0°); f\ - 1.256 GHz. 
Figure 3. Concluded. 
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(a) (Bi = 0°, (f>i = 90°); /o = 1.249 GHz (TM 0 i mode case). 

Figure 5. Vectorial representation of surface-current density on equilateral triangular patch (shown 
in fig. 2) excited by plane wave with angles of incidence of (0, =0°, 4>i = 90°) and W = 0 > 
6, = 0°), /o = 1.249 GHz (TM 0 i mode case), and /i = 1.276 GHz (TMio mode case). 




(b) (9i = 0°, <f>i = 0°); /i = 1.276 GHz (TM 10 mode case). 
Figure 5. Concluded. 
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Resonant frequency, GHz 



Resonant frequency, GHz 

Figure 6. Normalized current density as function of resonant frequency for equilateral triangular 
patch (shown in fig. 2) excited by plane wave at angle of incidence of (9 Z — 0° } fa = 0°) and 
M = N — 19. The terms f\ (1.276 GHz), fo (2.510 GHz), and /g (3.355 GHz) are resonant 
frequencies of TMio, TM 20 , and TM 21 modes, respectively. 
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Resonant frequency, GHz 


Figure 7. Normalized current density as function of resonant frequency for equilateral triangular 
patch (shown in fig. 2) excited by plane wave at angle of incidence of (0* = 0°, 4>i = 90°) and 
M = N = 19. The terms / 0 (1.249 GHz), f 2 (2.172 GHz), / 4 (2.525 GHz), and / 5 (3.265 GHz) 
are resonant frequencies of TM 0 i, TMu, TM 0 2, and TMi 2 modes, respectively. 
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(a) /2 = 2.172 GHz (TMjj mode case). 

Figure 8. Vectorial representation of surface-current density on equilateral triangular patch (shown in 
fig. 2) excited by plane wave with angle of incidence of (0* = 0°, — 0°) and resonant frequencies 

/2 — 2.172 GHz (TMn mode case), / 4 = 2.510 GHz (TM20 mode case), and = 3.355 GHz 
(TM21 mode case). 
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(b) /4 = 2.510 GHz (TM 20 mode case). 
Figure 8. Continued. 
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Figure 8. Concluded. 
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(a) h = 2.172 GHz (TMn mode case). 

Figure 9. Vectorial representation of surface-current density on equilateral triangular patch (shown in 
fig. 2) excited by plane wave with angle of incidence of {6 l = 0°, — 90°) and resonant frequencies 

fl — 2.172 GHz (TMn mode case), — 2.510 GHz (TMq 2 mode case), and = 3.265 GHz 
(TM 12 mode case). 
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(b) /4 = 2.510 GHz (TMq2 mode case). 


Figure 9. Continued. 
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(c) /g = 3.265 GHz (TM 12 mode case). 
Figure 9. Concluded. 
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12345ft## • • • • M + 1 

Figure 10. Geometry of circular patch with radius a = 1.88 cm, d = z' == 0.16 cm, e T — 2.53, and 
Loss tangent = 0.002. 
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(a) & = 0 °, & = 0°). 

Figure 11. Normalized current density as function of resonant frequency for circular patch (shown 
in fig. 10) for angles of incidence of (0j = 0°, fa — 0°) and (0j = 0°, fa = 90°), M = N = 11, and 
/o = 2.76 GHz. 
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(b) ($i = 0°, 4>i = 90°). 
Figure 11. Concluded. 
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11 12 13 14 15 16 17 18 19 20 2 1 22 23 

M = N 


Figure 12. Resonant frequency of circular patch TMq mode (shown in fig. 10), excited by plane 
wave at angle of incidence of (6 t = 0°, = 0°) as function of M = N. 


29 






Figure 13. Vectorial representation of surface-current density on circular patch (shown in fig. 10) 
excited by plane wave with angles of incidence of (0j = 0°, <fii = 0°) and (9 X = 0°, <j>i = 90°) and 
/o = 2.796 GHz (TMp mode case). 
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(b) (< 9 i - 0°, <t>i = 90°). 
Figure 13. Concluded. 
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(a) fi = 4.685 GHz (TM 21 mode case). 

Figure 14. Vectorial representation of surface-current density on circular patch (shown in fig. 10) 
excited by plane wave with angle of incidence of (0* = 45°, <j>i = 90°) and resonant frequencies 
/l = 4.685 GHz (TM 21 mode case), /-j = 5.855 GHz (TM 0 2 mode case), and / 3 = 6.360 GHz 
(TM 31 mode case). 


32 






iViV 




t 

I 

t 


t 

f 


/ : r 



r\~y Y~A~ 
'A X 


' V 


^ X*X'1 T 

k'k'k'K'A 
kxKk 

K- 


f 



V 


i 


f; 

A 


V: t 


•* 




T 











t- * 1 - 

L \ 








Y--Y-4-H- 


\ 


v 

^t 


t 


_\j._X_X>* 

VL^L.\f-Xi-X: 

* ' \i-WM-* 
U-U'-l 


V 


-iJ 


i 


(b) /2 = 5.855 GHz (TM 02 mode case). 
Figure 14. Continued. 



( c ) /3 = 6.360 GHz (TM 31 mode case). 


Figure 14. Concluded. 
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Figure 15. 
a = 45, 


a 

Geometry of trapezoidal patch with radius a = 8.08 cm, b = c = 4.04 cm, d = z f = 0.16 cm, 
e r = 2.65, and Loss tangent = 0.002. 






(a) & = 45°, 4>i = 0°). 

Figure 16. Normalized current density as function of resonant frequency for trapezoidal patch (shown 
in fig. 15) excited by plane wave at angles of incidence of ( 9 t = 45°, <f>i — 0°) and (9 t = 45° , 
— 90°), M = 20, and N - 10. 
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Figure 16. Concluded. 
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Figure 17. Vectorial representation of surface-current density on trapezoidal patch (shown in fig. 15) 
excited by plane wave with angle of incidence of (9 r = 0°, <f>i — 0°), /o = 1.345 GHz, M = 20, and 
N = 10. 
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